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In  phosphoric  acid  fuel  cell  (PAFC)  stack,  reduction  of  the  phosphoric  acid  that  is  impregnated  in  the 
cell  is  a  major  factor  to  restrict  the  operating  life.  In  this  paper,  the  phosphoric  acid  reduction  by  the 
evaporative  and  condensational  dissipation  on  the  effect  of  liquid  phase  migration  was  evaluated  by 
numerical  analysis.  The  calculations  that  include  the  behavior  of  evaporation,  condensation,  and  liquid 
phase  migration  for  the  phosphoric  acid,  and  the  behavior  of  gas  flow  and  heat  transfer  were  conducted  for 
the  model  cell  with  conditions  that  correspond  to  an  on  site  stack.  The  diffusion  coefficient  of  liquid  phase 
migration  was  derived  from  the  activation  energy  and  frequency  factor  based  on  the  experimental  results. 
The  phosphoric  acid  distributions  in  the  vapor  phase  and  liquid  phase  at  the  process  gas  and  electrode 
were  calculated  for  the  duration  of  3 136  h,  and  as  a  result  of  those,  the  exhausting  and  remaining  amount 
of  phosphoric  acid  for  the  cell  were  evaluated.  The  effects  of  liquid  phase  migration  on  the  phosphoric 
acid  dissipation  were  evaluated  by  analysis  for  various  diffusion  coefficients,  and  the  analysis  results 
were  compared  with  the  experimental  results  for  the  model  cell. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  phosphoric  acid  fuel  cell  (PAFC)  stack,  phosphoric  acid  that 
has  low  vapor  pressure  is  used  as  an  electrolyte  to  suppress  the 
evaporation  and  dissipation  with  the  process  gases  (the  fuel  and 
oxidant  gases).  However,  the  decrease  in  cell  performance  caused 
by  the  reduction  of  phosphoric  acid  is  not  negligible  at  the  oper¬ 
ating  temperature  of  about  200  °C.  For  stacks  that  are  designed  to 
have  an  operating  life  of  over  40,000  h,  a  basic  design  to  keep  low 
temperature  at  the  gas  exit  comparing  with  the  gas  inlet  of  the  cell 
is  adopted.  As  a  result  of  the  design,  the  evaporated  phosphoric  acid 
in  the  process  gas  is  condensed  and  recovered  to  the  cell  at  the  gas 
exit. 

Yoshioka  et  al.  have  conducted  experimental  and  numerical 
studies  for  evaporation  [1  ]  and  condensation  [2]  of  phosphoric  acid. 
In  those  studies,  the  evaporation  rate  was  treated  to  be  controlled 
by  the  diffusion  rate  in  the  concentration  boundary  layer  of  process 
gas  and  the  condensation  rate  was  treated  to  be  controlled  by  the 
cooling  rate  of  process  gas,  and  the  phosphoric  acid  distributions 
in  vapor  and  liquid  phases  and  the  recovery  of  phosphoric  acid  for 
the  model  cell  were  discussed. 
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In  the  previous  paper  [3],  numerical  analyses  that  include  the 
behavior  of  evaporation  and  condensation  for  the  phosphoric  acid, 
and  the  behavior  of  gas  flow  and  heat  transfer  were  conducted  for 
the  model  cell.  Time  series  solutions  for  the  phosphoric  acid  distri¬ 
butions  in  the  vapor  phase  and  liquid  phase  at  the  process  gas  and 
cell  were  calculated  duration  of  3300  h,  and  as  a  result  of  those, 
the  exhausting  and  remaining  phosphoric  acids  for  the  cell  were 
evaluated. 

In  this  paper,  the  liquid  phase  migration  of  phosphoric  acid  that 
has  been  eliminated  in  the  previous  paper  was  considered  with  the 
aim  of  investigation  for  more  practical  design  method  to  estimate 
the  precise  dissipation  of  phosphoric  acid.  The  diffusion  coefficient 
of  liquid  phase  migration  was  derived  from  the  activation  energy 
and  frequency  factor  based  on  the  experimental  results  for  the  sub¬ 
strate.  Numerical  analyses  that  include  the  behavior  of  evaporation, 
condensation,  and  liquid  phase  migration  for  the  phosphoric  acid, 
and  the  behavior  of  gas  flow  and  heat  transfer  were  conducted 
for  the  model  cell  that  consists  of  porous  electrode  substrate  and 
matrix  with  conditions  that  correspond  to  an  on  site  stack.  The 
phosphoric  acid  distributions  in  the  vapor  phase  and  liquid  phase 
and  the  exhausting  and  remaining  phosphoric  acids  for  the  model 
cell  were  calculated  duration  of  3 1 36  h,  and  as  a  result  of  those,  the 
exhausting  and  remaining  phosphoric  acids  for  the  cell  were  calcu¬ 
lated.  The  effects  of  liquid  phase  migration  on  the  phosphoric  acid 
dissipation  were  evaluated  by  calculations  for  the  various  diffusion 
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Nomenclature 

Ca 

phosphoric  acid  volume  (cm3) 

Cs 

substrate  pore  volume  (cm3) 

c 

phosphoric  acid  specific  volume 

D 

diffusion  coefficient  of  liquid  phase 
(cm2  s-1 ) 

migration 

Co 

initial  phosphoric  acid  specific  volume 

t 

time  (s) 

r 

radial  coordinate  (cm) 

Io 

modified  Bessel  function  of  1st  kind 

a 

phosphoric  acid  initial  dispersion  radius  1 

(cm) 

c  Ir=0 

central  phosphoric  acid  specific  volume 

J 

evaporated  or  condensed  phosphoric 
(mols-1  cm-2) 

acid  flux 

AG* 

activation  energy  change  (J  mol-1  K-1 ) 

c 

frequency  factor  for  J 

CF 

frequency  factor  for  diffusion  coefficient 

R 

universal  gas  constant  (8.314J  mol-1 1<-1 

) 

T 

absolute  temperature  (K) 
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Fig.  1.  Outline  of  liquid  phase  migration  experiment. 


coefficients,  and  the  analysis  results  were  evaluated  by  comparing 
with  the  experimental  results. 

2.  Experimental  methods  and  results 

2. 2.  Liquid  phase  migration 


1  2  3  4  5 

Sample  Piece  No. 

Fig.  2.  Specific  volume  distributions  after  1  day. 


the  substrate  pore  volume  Cs,  the  phosphoric  acid  specific  volume 
C  was  defined  as  follows. 

c=£  (1) 

Us 

Regarding  to  the  time  variable  specific  volume  diffusion  for  the 
cylindrical  symmetric  one-dimensional  direction  that  represents 
the  radial  coordinate  of  the  substrate,  the  specific  volume  distribu¬ 
tion  was  derived  as  follows. 


Here,  Co  is  initial  specific  volume  of  phosphoric  acid  at  the  sub¬ 
strate  center  and  equals  to  1.0  in  each  experiment,  D  is  diffusion 
coefficient  cm2  s-1,  t  is  time  s,  r  is  radial  coordinate  cm,  r'  is  radial 
coordinate  cm,  /0  is  modified  Bessel  function  of  1st  kind,  and  a  is 
initial  dispersion  radius  cm  of  phosphoric  acid. 

The  specific  volume  ratio  distributions  i.e.  the  specific  volume 
distributions  (C/C0  =  C)  of  the  sample  pieces  for  various  tempera¬ 
tures  are  shown  in  Fig.  2.  Based  on  these  results,  the  central  specific 
volumes  C|r=0  were  derived  by  extrapolation  for  each  case. 

Regarding  to  the  specific  volume  variation  at  the  center  of  the 
cylindrical  coordinate,  the  following  equation  was  derived  from  Eq. 
(2)  by  substituting  0.0  for  r. 


The  outline  of  the  migration  experiment  is  shown  in  Fig.  1,  and 
the  experimental  conditions  are  shown  in  Table  1.  The  substrate 
size  was  20  cm  x  20  cm  each,  and  the  phosphoric  acid  that  occupied 
10%  of  the  substrate  pore  volume  (3.67  cc  of  85wt%  H3P04,  25  °C) 
was  impregnated  and  dispersed  uniformly  by  keeping  in  50  °C  air 
for  5  days.  After  that,  the  phosphoric  acid  (1.0  cc  of  85  wt%  H3PO4, 
25  °C)  was  dropped  at  the  center  of  substrate  using  an  injection 
syringe  and  the  initial  dispersion  radius  was  measured  for  each 
substrate.  Each  substrate  was  kept  in  various  temperature  con¬ 
ditions  as  shown  in  Table  1  for  1  day,  and  the  phosphoric  acid 
concentrations  in  each  sample  piece  shown  in  Fig.  1  were  mea¬ 
sured  using  ICP-OES  (Inductively  Coupled  Plasma  Optical  Emission 
Spectrometer,  ICPS-5000/Shimadzu).  The  phosphoric  acid  concen¬ 
trations  were  converted  to  the  phosphoric  acid  volume  (85wt% 
H3P04,  25  °C),  and  based  on  the  phosphoric  acid  volume  Ca  and 

Table  1 

Conditions  for  liquid  phase  migration  experiments. 

Temperature 
Pressure 
Surrounding 
Duration  time 
Fallen  phosphoric  acid 
Initial  dispersion  radius 


c  =  c°{'-exp(ro)}  (3) 

The  diffusion  coefficients  D  were  derived  by  Eq.  (3)  according  to 
C|r=0,  a,  and  t  for  each  temperature  condition.  The  derived  diffusion 
coefficients  are  shown  in  Fig.  3,  and  the  following  equation  was 
derived  by  regression  of  these  results  to  the  Arrhenius  formula. 


2.2.  Evaporative  and  condensational  dissipation 

The  outline  of  the  model  cell  dissipation  experiment  is  shown 
in  Fig.  4,  and  the  experimental  conditions  are  shown  in  Table  2.  The 
experiment  was  carried  out  using  a  porous  electrode  substrate  and 
matrix  in  that  the  phosphoric  acid  was  impregnated,  and  the  steam 
humidified  N2  gas  that  substitute  process  gas  was  run  through  the 
electrode  gas  channel.  The  temperature  in  each  electrode  section 
was  controlled  by  electrical  heaters  to  correspond  with  the  cell 
temperature  distribution  along  the  oxidant  gas  flow  for  an  on  site 
stack.  The  liquid  phase  phosphoric  acid  was  able  to  migrate  through 
the  electrode  sections,  as  there  were  no  barriers  between  the  elec¬ 
trode  sections.  The  amount  of  the  remaining  phosphoric  acid  in 


50, 100, 165,  and  205  °C 
1  atm 

Air  +  60  °C  humidified 
1  day 

1  cc  (85  wt%  H3PO4,  25  °C) 
2. 1-2.8  cm 


D  =  6.89  x  10  3  x  exp 


-1.50  x 
RT 
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Fig.  3.  Diffusion  coefficient  of  liquid  phase  migration. 
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Fig.  4.  Outline  of  model  cell  dissipation  experiment. 


Table  2 

Conditions  for  dissipation  experiment  and  analysis. 

Pressure 

Gas  flow  rate 

Operation  duration 

Initial  phosphoric  acid  amount 

1  atm 

N2  (8.3 1  min-1 )  +  60  °C  humidified 
3136h 

9.6 cc  (85  wt%  H3PO4,  25  °C)/section 

Temperature  Inlet  Section 

Outlet 

1  2 

3  4  5 

205  °C  205  °C 

201 °C  179°C  163  °C 

(b)  Migration  in  the  liquid  phase.  That  is  migration  by  capillary 
force  caused  by  the  concentration  or  temperature  change  of 
phosphoric  acid  in  the  porous  electrode. 


In  the  previous  paper,  the  phenomenon  above  (a)  was  consid¬ 
ered  on  the  calculation.  The  phosphoric  acid  was  regarded  as  a  two 
components  system  of  phosphorus  pentoxide  (P4O10)  and  water 
(H20),  and  the  evaporation  and  condensation  rates  of  phosphoric 
acid  were  derived  by  the  following  activation  process  based  on  the 
vapor  liquid  equilibrium  of  a  two  components  system  in  reference 
[5-8]  and  the  nucleation  theory  of  a  two  components  system  in 
reference  [9-12]. 


J  =  c  exp 


(5) 


Here,  J  is  the  amount  of  evaporating  or  condensing  phospho¬ 
ric  acid  per  unit  area  and  unit  time  on  the  electrode  surface 
mols-1  cm-2,  R  is  the  universal  gas  constant  8.314Jmol-1  I<-1,  T 
is  the  liquid  phase  absolute  temperature  K,  AG*  is  the  activation 
energy  J  mol-1,  and  c  is  the  frequency  factor. 

In  this  paper,  besides  of  the  above  phenomenon  (a),  the  phe¬ 
nomenon  (b)  was  considered,  and  the  balance  of  phosphoric  acid 
of  liquid  phase  in  the  electrode  was  expressed  as  follows. 


dC  ( d2c  d2c\ 
~dt  ~  Ivfo2  +  sy2y/ 


(6) 


Here,  t  is  time,  x  and  y  are  coordinates  on  the  electrode  surface,  C  is 
specific  volume  of  the  phosphoric  acid,  and  D  is  diffusion  coefficient 
derived  by  Eq.  (4)  as  a  function  of  T. 

The  unsteady  analysis  loop  that  includes  the  analysis  loop  of 
the  code  for  MCFC  was  programmed.  The  behaviors  of  evaporation, 
condensation,  liquid  phase  migration  of  the  phosphoric  acid  were 
calculated,  and  as  results  of  these,  the  amount  of  phosphoric  acid  in 
the  vapor  phase  and  liquid  phase  at  the  process  gas  and  electrode 
were  calculated  in  the  unsteady  loop. 

Analyses  were  conducted  on  the  conditions  that  correspond  to 
the  model  cell  experiment  that  was  mentioned  in  Section  2.2,  using 
a  one  dimensional  model  in  the  direction  of  the  gas  flow  and  each 
electrode  section  was  divided  into  nine  elements. 


4.  Analysis  results 


each  electrode  section  was  derived  from  the  quantitative  analysis 
results  by  ICP-OES  after  3136  h  elapsed  state, 

3.  Analysis  method  and  conditions 

For  the  calculations,  the  analysis  code  that  was  developed  for  a 
PAFC  stack  [3]  was  used  with  the  addition  of  a  function  for  the  liq¬ 
uid  phase  migration.  The  previous  code  was  based  on  the  analysis 
code  that  has  been  developed  for  the  characteristic  calculation  of  a 
molten  carbonate  fuel  cell  (MCFC)  stack  [4],  and  the  function  of  elec¬ 
trode  reaction  has  modified  to  correspond  to  that  of  the  PAFC  and 
the  functions  for  the  evaporation  and  condensation  of  phosphoric 
acid  have  been  added. 

The  phosphoric  acid  migrations  in  a  PAFC  stack  are  considered 
to  be  caused  by  the  following  two  mechanisms. 

(a)  Migration  accompanied  with  the  phase  change  between  liq¬ 
uid  and  vapor.  That  is  dissipation  of  phosphoric  acid  from  the 
porous  electrode  to  process  gas  by  evaporation  and  recovery 
from  the  process  gas  to  porous  electrode  by  condensation  with 
the  transfer  by  gas  flow. 


4.1.  Amount  of  phosphoric  acid  in  liquid  phase  at  the  electrode 

The  amount  of  remaining  P4O10  in  each  electrode  section  at 
3 1 36  h  elapsed  state  for  the  analysis  results  and  experimental  result 


Electrode  Section  No. 
Fig.  5.  Residual  P4O10  in  the  electrode. 
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are  shown  in  Fig.  5.  These  were  derived  from  the  evaporation  and 
condensation  of  P4O10  and  H20  and  the  liquid  phase  migration 
of  phosphoric  acid  that  will  be  mentioned  later.  Here,  besides  of 
the  diffusion  coefficient  shown  in  Eq.  (4),  the  following  diffusion 
coefficients  were  adopted  for  the  analysis. 


0 

II 

0 

(7) 

D  =  2.00  x  10  2  x  exp  | 

—1.50  x  104\ 

(8) 

.  RT  ) 

D  =  3.00  x  10-2  x  exp  | 

—1.50  x  104\ 

(9) 

RT  ) 

As  shown  in  Fig.  5,  there  was  little  agreement  between  the  anal¬ 
ysis  result  derived  from  the  diffusion  coefficient  of  Eq.  (4)  and  the 
experimental  result,  though  the  analysis  result  derived  from  the 
diffusion  coefficient  of  Eq.  (9)  was  in  good  agreement  with  the 
experimental  result.  Here,  it  should  be  noted  that  the  electrode 
substrate  and  matrix  were  used  for  the  dissipation  experiments; 
on  the  other  hand,  the  electrode  substrate  was  used  for  the  liquid 
phase  migration  experiments.  As  mentioned  in  [13],  typical  sub¬ 
strate  is  around  370  |xm  thickness  and  around  80%  porosity,  though 
typical  matrix  is  around  50  [xm  thickness  and  lesser  porosity.  There¬ 
fore,  in  the  electrode,  the  phosphoric  acid  is  mainly  reserved  in 
the  matrix  that  has  lesser  porosity  comparing  with  the  substrate 
and  the  liquid  phase  migration  characteristic  was  considered  to  be 
mainly  controlled  by  the  matrix.  In  other  words,  the  diffusion  coef¬ 
ficient  according  to  the  electrode  consists  of  substrate  and  matrix 
was  considered  to  be  different  from  the  diffusion  coefficient  accord¬ 
ing  to  the  substrate.  Namely,  it  was  considered  that  the  diffusion 
coefficient  of  Eq.  (4)  showed  that  of  the  substrate,  and  the  diffu¬ 
sion  coefficient  of  Eq.  (9)  showed  that  of  the  electrode  consists  of 
substrate  and  matrix. 

Consequently,  the  analysis  result  according  to  the  diffusion 
coefficient  of  Eq.  (9)  was  adopted  for  the  comparison  with  the 
experimental  result.  The  analysis  results  according  to  the  diffusion 
coefficients  of  Eqs.  (7)  and  (9)  were  discussed  regarding  the  effect 
of  the  liquid  phase  migration  that  will  have  an  influence  on  the 
phosphoric  acid  dissipation.  Here,  the  analysis  results  according  to 
the  diffusion  coefficient  of  Eq.  (7)  show  the  results  eliminating  the 
liquid  phase  migration. 

In  Fig.  5,  after  3136  h,  the  amount  of  remaining  P4O10  for  the 
analysis  and  experimental  results  increased  according  to  the  sec¬ 
tion  goes  to  the  gas  outlet  direction  in  each  case.  In  sections  1 
and  2,  P4O10  decreased  from  the  impregnated  P4O10  at  the  initial 
state,  and  the  dissipation  of  phosphoric  acid  to  the  process  gas  by 
evaporation  was  observed.  In  sections  4  and  5,  there  were  more 
P4O10  than  the  impregnated  P4O10,  and  the  recovery  of  phospho¬ 
ric  acid  from  the  process  gas  by  condensation  was  observed.  The 
decreases  of  remaining  P4O10  in  sections  1  and  2,  and  the  increases 
in  section  4  and  5  became  lesser  in  accordance  with  the  increase 
of  liquid  phase  migration  diffusion  coefficient.  These  showed  that 
the  liquid  phase  migration  deregulated  the  dissipation  and  recov¬ 
ery  of  P4O10  by  the  evaporation  and  condensation.  As  a  result  of 
those,  the  residual  P4O10  decreased  0.00836  mol  in  section  1,  and 
increased  0.0061 6  mol  in  section  5  for  the  diffusion  coefficient  of 
Eq.  (9),  though  the  residual  P4O10  decreased  0.0347  mol  in  section 
1  and  increased  0.0620  mol  in  section  5  for  the  diffusion  coefficient 
of  Eq.  (7). 

The  total  residual  P4O10  in  the  electrode  for  the  various  diffusion 
coefficient  frequency  factors  cF  are  shown  in  Fig.  6.  The  exhausted 
P4O10  after  3136  h  were  0.00475  mol  and  0.00558  mol  for  the  dif¬ 
fusion  coefficient  of  Eqs.  (7)  and  (9),  respectively.  As  shown  here, 
the  difference  of  total  exhausted  P4O10  for  various  liquid  phase 
migration  diffusion  coefficients  were  a  few.  On  the  other  hands, 
as  mentioned  about  Fig.  5,  the  difference  of  decreasing  P4O10  in 


Frequency  Factor  of  Diffusion  Coefficient 
cF 

Fig.  6.  Total  residual  P4O10  in  the  electrode. 


0  0.2  0.4  0.6  0.8  1 

Distance  from  Inlet  (m) 

Fig.  7.  P4O10  and  H20  molar  density  of  liquid  phase.  In  the  electrode,  D  =  0.0. 


sections  1  and  2  for  various  liquid  phase  migration  diffusion  coeffi¬ 
cients  were  much  and  the  cell  performance  decline  in  these  portion 
was  considered  to  be  fatal  in  the  case  of  eliminated  liquid  phase 
migration.  Consequently,  the  liquid  phase  migration  had  a  big  effect 
on  the  P4O10  dissipation  and  on  the  cell  performance  through  P4O10 
decrease  in  the  sections  of  gas  inlet  side. 

P4O10  and  H20  molar  concentration  distributions  in  liquid  phase 
at  electrode  along  the  gas  flow  direction  at  the  initial,  2000  h,  and 
3136  h  elapsed  state  for  the  diffusion  coefficients  of  Eqs.  (7)  and  (9) 
are  shown  in  Figs.  7  and  8,  respectively. 
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Fig.  8.  P4O10  and  H20  molar  density  of  liquid  phase.  In  the  electrode, 
D  =  3.00 E  -  2  exp( - 1 .50F4/KT). 
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Fig.  9.  Evaporation  and  condensation  rates  along  the  electrode,  minus:  evaporation, 
plus:  condensation,  D  =  0.0. 
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Fig.  10.  Evaporation  and  condensation  rates  along  the  electrode,  minus:  evapora¬ 
tion,  plus:  condensation,  D  =  3.00E- 2  exp(-1.50E4/i?r). 


In  the  case  for  diffusion  coefficient  of  Eq.  (7),  the  molar  con¬ 
centration  of  P4O10  and  H20  decreased  at  sections  1  and  2  these 
are  0.0-0.4  m  distance  from  the  gas  inlet,  and  increased  at  section 
5  that  is  0.8-1. 0  m  distance  from  the  gas  inlet  with  the  progress 
of  time.  Here,  the  decrease  and  increase  of  P4O10  and  H20  molar 
concentration  were  caused  by  the  evaporation  and  condensa¬ 
tion  accompanying  gas  flow  transform,  respectively.  Since  P4O10 
molar  fraction  of  evaporating  and  condensing  phosphoric  acid  were 
thicker  than  that  of  impregnated  phosphoric  acid  (0.0922;  85  wt% 
H3PO4)  in  the  electrode  as  it  will  be  shown  in  Fig.  11,  the  varia¬ 
tion  of  P4O10  concentration  in  liquid  phase  at  electrode  was  larger 
than  the  variation  of  H20  concentration  at  sections  1,  2,  and  5.  It 
will  be  shown  in  Fig.  9,  the  peak  position  of  evaporation  and  con¬ 
densation  rates  moved  toward  the  gas  downstream  and  these  rates 
decreased  with  the  time  progress.  These  were  correspondent  with 
the  decrease  and  increase  of  the  molar  concentration  of  P4Oi0  and 
H20  in  the  liquid  phase  according  to  the  gas  flow  for  the  evaporation 
and  condensation,  respectively,  in  Fig.  7. 

In  the  case  for  diffusion  coefficient  of  Eq.  (9),  the  molar  con¬ 
centration  of  P4O10  and  H20  decreased  at  sections  1  and  2,  and 
increases  at  section  5  with  the  progress  of  time.  In  this  case,  the 
change  for  molar  concentration  of  P4O10  and  H20  in  each  section 
was  lesser  than  the  case  for  the  diffusion  coefficient  of  Eq.  (7).  These 
were  caused  by  the  liquid  phase  phosphoric  acid  recirculation  from 
section  5  to  sections  1  and  2  by  liquid  phase  migration,  and  as  result 
of  that,  the  liquid  phase  conditions  were  kept  almost  stable  through 
the  3136  h  duration.  On  the  other  hands,  the  liquid  phase  phospho¬ 
ric  acid  conditions  were  continuously  changed  with  the  progress  of 
time,  as  the  liquid  phase  migration  was  eliminated  in  the  case  of 
Fig.  7. 

4.2.  Evaporation  and  condensation  rates  of  phosphoric  acid 

The  evaporation  and  condensation  rate  distributions  of  phos¬ 
phoric  acid  along  the  gas  flow  direction  for  the  initial,  2000  h,  and 
3136  h  elapsed  state  for  the  diffusion  coefficient  of  Eqs.  (7)  and  (9) 
are  shown  in  Figs.  9  and  10,  respectively.  Here,  the  negative  rate 
shows  evaporation  and  positive  rate  shows  condensation  in  each 
figure. 

In  the  case  for  diffusion  coefficient  of  Eq.  (7),  the  evaporation 
arose  at  sections  1  and  2,  and  the  condensation  arose  at  section  5 
where  the  electrode  temperature  steeply  fell  comparing  with  the 
upstream.  The  peak  magnitude  of  evaporation  and  condensation 
rate  fell  in  accordance  with  the  progress  of  time,  and  the  peak  posi¬ 
tion  of  evaporation  and  condensation  rate  gradually  moved  to  the 
downstream  of  gas  flow. 

In  the  case  for  diffusion  coefficient  of  Eq.  (9),  at  the  initial  state, 
the  evaporation  and  condensation  rate  distribution  was  the  same 
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Fig.  11.  P4O10  molar  fractions  of  evaporating  and  condensing  phosphoric  acid  along 
the  electrode,  D  =  0.0. 


as  the  case  for  the  diffusion  coefficient  of  Eq.  (7).  However,  the  peak 
magnitude  and  position  of  evaporation  and  condensation  rate  were 
almost  stable  in  the  progress  of  time.  These  showed  that  the  liquid 
phase  phosphoric  acid  migrated  from  section  5  to  sections  1  and  2 
in  the  progress  of  time,  and  were  correspondent  with  the  difference 
between  the  results  of  Figs.  7  and  8. 

P4O10  molar  fraction  distribution  of  evaporating  and  condensing 
phosphoric  acid  composition  along  the  gas  flow  direction  for  the 
diffusion  coefficients  of  Eqs.  (7)  and  (9)  are  shown  in  Figs.  1 1  and  1 2, 
respectively. 
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Fig.  12.  P4O10  molar  fractions  of  evaporating  and  condensing  phosphoric  acid  along 
the  electrode,  D  =  3.00E  -  2  exp( - 1 .50E4/R7). 
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Fig.  13.  P4O10  partial  pressure  along  the  electrode,  D  =  0.0. 

In  the  cases  for  diffusion  coefficient  of  Eq.  (7),  P4O10  molar 
fraction  of  evaporating  composition  at  sections  1,  2,  and  3  fell 
in  accordance  with  the  time  progress,  and  P4O10  molar  fraction 
changed  from  0.745  to  0.743  in  section  1  through  the  3136  h  dura¬ 
tion.  In  condensing  composition  at  section  5,  P4O10  molar  fraction 
increased  in  according  with  the  time  progress,  and  changed  from 
0.747  to  0.752  through  the  3136  h  duration. 

In  the  cases  for  diffusion  coefficient  of  Eq.  (9),  P4O10  molar  frac¬ 
tion  change  of  evaporating  and  condensing  compositions  in  the 
time  progress  were  lesser  than  the  cases  for  diffusion  coefficient 
of  Eq.  (7),  and  P4O10  molar  fraction  change  was  from  0.7454  to 
0.7450  at  section  1  and  from  0.7454  to  0.7457  at  section  5  through 
the  3136  h  duration. 

P4O10  concentration  of  evaporating  and  condensing  composi¬ 
tions  were  thicker  than  the  molar  fraction  0.0922  (85  wt%  H3PO4, 
25  °C)  of  impregnated  phosphoric  acid  for  each  case. 

4.3.  Partial  pressure  and  saturated  vapor  pressure  of  phosphoric 
acid 

P4O10  partial  pressure  distributions  of  vapor  phase  in  the  pro¬ 
cess  gas  along  the  gas  flow  direction  at  the  initial,  2000  h  and  3 1 36  h 
elapsed  states  for  the  diffusion  coefficients  of  Eqs.  (7)  and  (9)  are 
shown  in  Figs.  13  and  14,  respectively. 

In  the  case  for  diffusion  coefficient  of  Eq.  (7),  at  the  initial 
state,  P4O10  partial  pressure  increased  gradually  along  the  gas  flow 
in  sections  1  and  2  where  evaporation  occurred,  and  decreased 
steeply  and  reached  to  fixed  values  in  section  5  where  condensation 
occurred.  The  gradual  increases  at  section  1  and  steep  decrease  at 
section  5  for  P4O10  partial  pressures  were  correspondent  with  the 
gradual  decrease  of  the  evaporation  rate  at  section  1  and  the  steep 
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Fig.  14.  P4O10  partial  pressure  along  the  electrode,  D  =  3.00£-  2  exp(-1.50£4/K'T). 


Fig.  15.  P4O10  saturated  vapor  pressure  and  partial  pressure  along  the  electrode, 
D  =  0.0. 

increase  of  condensation  rate  at  section  5  in  Fig.  9.  The  increase  of 
P4O10  partial  pressures  at  the  gas  inlet  part  became  lesser  accord¬ 
ing  to  the  progress  of  time,  and  the  distance  between  the  gas  inlet 
and  the  position  where  the  vapor  pressure  reached  to  saturation 
became  long.  P4O10  partial  pressure  at  the  gas  exit  that  determines 
the  amount  of  exhausting  phosphoric  acid  slightly  increased  from 
0.00603  Pa  to  0.00654  Pa  in  the  duration  of  3 1 36  h. 

In  the  case  for  diffusion  coefficient  of  Eq.  (9),  at  the  initial  state, 
P4O10  partial  pressure  distribution  was  the  same  as  the  case  for  the 
diffusion  coefficient  of  Eq.  (7).  The  changes  of  P4Oi0  partial  pres¬ 
sures  distributions  according  to  the  progress  of  time  were  lesser 
than  the  case  for  diffusion  coefficient  of  Eq.  (7)  and  P4O10  partial 
pressure  distributions  were  almost  stable  except  sections  3  and  4  in 
the  duration  of  3136  h.  These  showed  the  effect  of  the  phosphoric 
acid  liquid  phase  migration  from  section  5  to  sections  1  and  2,  and 
were  correspondent  with  the  steadiness  for  the  evaporation  and 
condensation  rates  with  the  progress  of  time  as  shown  in  Fig.  10. 
P4O10  partial  pressure  was  around  0.00603  Pa  on  the  gas  exit  at 
3136  h  elapsed  state,  and  it  was  almost  same  as  that  of  the  case  for 
diffusion  coefficient  of  Eq.  (7).  This  caused  the  total  residual  P4O10 
to  be  almost  same  between  the  cases  for  diffusion  coefficients  of 
Eqs.  (7)  and  (9)  as  shown  in  Fig.  6. 

P4O10  partial  pressure  in  the  vapor  phase  and  P4O10  saturated 
vapor  pressure  that  was  determined  by  the  condition  in  liquid 
phase  at  electrode  at  initial  state  and  3136  h  elapsed  state  for  the 
diffusion  coefficients  of  Eqs.  (7)  and  (9)  are  shown  in  Figs.  1 5  and  1 6, 
respectively. 

In  the  case  for  diffusion  coefficient  of  Eq.  (7),  at  the  initial  state, 
P4O10  partial  pressure  at  the  gas  inlet  part  was  lower  than  the 
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Fig.  16.  P4O10  saturated  vapor  pressure  and  partial  pressure  along  the  electrode, 
D  =  3.00 E  -  2  exp( - 1 .50£4/£r). 
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saturated  vapor  pressure,  and  according  to  the  gas  flow  the  par¬ 
tial  pressure  increased  and  reached  to  saturation.  At  sections  3  and 

4  these  are  0.4-0.8  m  distance  from  the  gas  inlet,  P4O10  partial  pres¬ 
sure  was  higher  than  the  saturated  vapor  pressure  and  the  vapor 
phase  of  process  gas  was  supersaturating  state  regarding  to  the 
temperature  and  phosphoric  acid  concentration  in  the  liquid  phase 
at  the  electrode.  At  section  5,  P4O10  partial  pressure  decreased 
steeply  and  reached  to  saturated  vapor  pressure  according  to  the 
liquid  phase  conditions  at  the  gas  exit. 

At  the  3136  h  elapsed  state,  P4O10  saturated  vapor  pressure 
increased  from  0  according  to  the  gas  flow,  and  P4O10  partial  pres¬ 
sure  distribution  was  almost  equal  to  the  saturated  vapor  pressure 
distribution  in  sections  1  and  2.  P4O10  saturated  vapor  pressure 
decreased  comparing  with  the  initial  state,  as  the  liquid  phase  phos¬ 
phoric  acid  concentration  was  decreased  because  of  the  thicker 
evaporating  P4O10  concentration  as  shown  in  Fig.  11.  P4O10  par¬ 
tial  pressure  was  supersaturating  state  at  section  4,  and  decreased 
gradually  at  section  5  to  reach  to  the  saturated  vapor  pressure  at 
the  gas  exit.  The  saturated  vapor  pressure  steeply  rose  at  the  gas 
inlet  portion  of  section  5  because  of  the  steep  increase  of  P4O10 
concentration  in  the  electrode  by  condensation  as  shown  in  Fig.  7. 

In  the  case  for  diffusion  coefficient  of  Eq.  (9),  at  the  initial  state, 
P4O10  partial  pressure  and  saturated  vapor  pressure  distributions 
were  the  same  as  the  case  for  the  diffusion  coefficient  of  Eq.  (7)  in 
each. 

At  the  3136  h  elapsed  state,  P4O10  saturated  vapor  pressure 
increased  from  0  according  to  the  gas  flow,  and  P4O10  partial  pres¬ 
sure  distribution  was  lower  than  the  saturated  vapor  pressure 
distribution  in  sections  1  and  2.  P4O10  saturated  vapor  pressure 
was  lower  comparing  to  the  initial  state,  however  the  difference 
was  lesser  than  the  case  for  diffusion  coefficient  of  Eq.  (7).  The  dif¬ 
ference  of  P4O10  partial  pressure  distributions  between  the  initial 
state  and  3136  h  elapsed  state  was  also  lesser  than  the  case  for 
diffusion  coefficient  of  Eq.  (7).  P4O10  partial  pressure  was  supersat¬ 
urating  state  at  sections  3  and  4,  and  it  decreased  steeply  at  section 

5  to  reach  to  the  saturated  vapor  pressure  at  the  gas  exit.  The  sat¬ 
urated  vapor  pressure  at  section  5  was  almost  same  through  the 
duration  of  3136  h;  this  was  caused  by  the  steadiness  of  phospho¬ 
ric  acid  in  section  5  by  the  liquid  phase  migration  from  section  5 
to  sections  1  and  2.  This  was  correspondent  with  the  steadiness 
in  the  progress  of  time  for  the  evaporation  and  condensation  rates 
and  liquid  phase  molar  densities  of  phosphoric  acid  as  shown  in 
Figs.  10  and  8,  respectively. 

5.  Conclusion 

The  dissipation  characteristics  of  phosphoric  acid  in  PAFC  were 
evaluated  by  the  analysis  using  a  code  that  includes  the  behavior 
of  the  evaporation,  condensation,  and  liquid  phase  migration  of 
phosphoric  acid  and  gas  flow  and  heat  transfer.  The  analysis  was 
conducted  for  the  model  cell  consists  of  porous  electrode  substrate 
and  matrix  that  has  correspondent  conditions  for  an  on  site  stack 
in  the  duration  of  3136  h. 


The  liquid  phase  migration  of  phosphoric  acid  were  calculated 
using  an  Arrhenius  formula  diffusion  coefficients  derived  from 
experimental  result  for  substrate,  and  the  phosphoric  acid  dissi¬ 
pation  characteristic  for  each  electrode  section  was  calculated.  The 
analysis  results  had  good  agreement  with  the  experimental  results 
by  adopting  appropriate  frequency  factor  for  diffusion  coefficient. 
Consequently  the  effectiveness  of  the  present  analysis  method  was 
validated. 

The  effects  of  liquid  phase  migration  on  the  phosphoric  acid 
dissipation  were  evaluated  by  comparing  the  analysis  results 
including  the  liquid  phase  migration  and  eliminating  that.  In  the 
case  of  eliminating  liquid  phase  migration,  the  amount  of  remain¬ 
ing  P4O10  decreased  in  the  high  temperature  sections  at  gas  inlet 
side  and  increased  in  the  low  temperature  sections  at  gas  outlet 
side  through  the  3136  h  duration.  Namely,  P4O10  in  the  gas  inlet 
sections  migrated  to  the  gas  outlet  sections  by  the  evaporation  and 
condensation  accompanied  with  gas  flow.  On  the  other  hand,  in 
the  case  of  including  liquid  phase  migration,  P4O10  decrease  and 
increase  in  the  high  and  low  temperature  sections  were  lesser  than 
the  above  case.  As  a  result  of  that,  P4O10  migrations  from  the  gas 
inlet  sections  to  the  gas  outlet  sections  were  lesser  than  the  case  of 
eliminating  liquid  phase  migration.  Consequently,  it  was  observed 
that  the  liquid  phase  migration  deregulated  the  effect  of  migra¬ 
tion  by  the  evaporation  and  condensation  accompanied  with  gas 
flow. 

Both  cases  of  including  and  eliminating  liquid  phase  migra¬ 
tion,  P4O10  partial  pressure  was  affected  by  the  conditions  of  the 
vapor  phase  in  process  gas  and  the  liquid  phase  in  electrode,  and 
for  some  conditions,  P4O10  partial  pressure  could  be  higher  than 
the  saturated  vapor  pressure  according  to  the  liquid  phase  in  the 
electrode,  and  the  supersaturating  states  were  resulted  to  be  pos¬ 
sible. 

As  a  result  of  these,  the  remaining  distribution  and  exhausting 
amount  of  P4Oi0  for  the  electrode  were  estimated  in  detail. 
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